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Abstract. The local purity of large many-body quantum systems can be studied 
by following a statistical mechanical approach based on a random matrix model. 
Restricting the analysis to the case of global pure states, this method proved to be 
successful and a full characterization of the statistical properties of the local purity was 
obtained by computing the partition function of the problem. Here we generalize these 
techniques to the case of global mixed states. In this context, by uniformly sampling 
the phase space of states with assigned global mixedness, we determine the exact 
expression of the first two moments of the local purity and a general expression for 
the moments of higher order. This generalizes previous results, obtained for globally 
pure configurations. Furthermore, through the introduction of a partition function 
for a suitable canonical ensemble, we compute the approximate expression of the first 
moment of the marginal purity in the high temperature regime. In the process, we 
establish a formal connection with the theory of quantum twirling maps that provides 
an alternative, possibly fruitful, way of performing the calculation. 
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1. Introduction 

In quantum mechanics the purity of a quantum state measures the amount of 
indeterminacy which originates either from lack of information on the process that 
led to its preparation, or from the correlations that have been established with some 
external party (environment). Mathematically speaking the two mechanisms are 
strongly interconnected, as the former can always be represented as an instance of 
the latter via purification [H [2]. For composite quantum system one can distinguish 
between global purity (i.e. the purity of the system as a whole), and the local (or 
"marginal") purities (i.e the purities associated with the subsystems that compose it). 
Clarifying the connection between the global and local purities of a quantum state 
is an important problem of quantum information theory [31 IH El E] which is closely 
related to the characterization of bipartite entanglement [3 IE]- In particular, for 
pure global states, the local purity of a subsystem provides a direct measure of the 
bipartite entanglement between the two parts: the smaller the purity, the larger the 
bipartite entanglement. For mixed (non pure) global states instead, the connection 
between local purity and entanglement is more subtle: no direct relation between the 
two quantities exists and bipartite entanglement measures for the global system can 
be obtained only by taking proper convex roof extensions of the local purity (the 
average being computed over all convex decompositions of the initial global state). Still, 
studying how the global indeterminacy of a composite system affects the indeterminacy 
of its constituents is important on its own and raises fundamental theoretical questions 
which are deeply interwoven with thermodynamical issues and call for a thoughtful 
investigation. Specifically here we focus on the following basic question: given a certain 
level of global mixedness of a many-body quantum system (say a mixture of gases at 
thermal equilibrium at a given temperature), which portion of such indeterminacy can be 
"accounted for" by its constituents? (that is how much the local mixedness contributes 
to the global one?). While for classical systems global and local mixedness are strongly 
interconnected (e.g., the former is always greater than the latter), in quantum mechanics 
the relation is more ambiguous, as qualitatively different sort of correlations can be 
established between the various subsystems. 

In order to account for the large spectrum of possibilities, we resort to a statistical 
approach, by exploiting tools and techniques imported from classical statistical 
mechanics. The distribution of the local purity of pure global states of large composite 
systems was studied in Refs. (91 |TQ] by using statistical mechanical methods. Here we 
endeavor to generalize the same techniques to a system in a global mixed state. In 
particular, by uniformly sampling the states with a given value of the global purity, we 
determine the exact expression of the first two moments of the local purity and obtain a 
general formula, valid for arbitrary moments. These findings generalize previous results 
obtained for globally pure configurations by Lubkin [15], Page [T7j, Lloyd and Pagels [16], 
Scott and Caves [32], and Giraud [2"HI2"5] . Then, through the introduction of a partition 
function for a suitable canonical ensemble endowed with a Lagrange multiplier, that 



Statistical distribution of the local purity in a large quantum system 



3 



plays the role of a fictitious temperature, we compute the approximate expression of 
the first moment of the marginal purity in the high temperature regime. We also 
establish the scaling of all these quantities with the dimension of the system, in the 
thermodynamic limit. Furthermore, an interesting connection appears between our 
problem and the theory of quantum channels. More precisely, the symmetry properties 
of the twirling transformations [26] can be proved to be very useful for the computation 
of the exact expression of the first moment of the local purity. As a final remark, let us 
stress the key role played by the introduction of a partition function for the system: it 
enables us to translate our problem, and thus the results we found, in terms of the even 
more general context of the theory of random matrices. 

This paper is organized as follows. In Sec. [2] we introduce the notation and set the 
basis of the statistical mechanical approach to the problem, starting from the simpler 
case of pure states and generalizing it to the case of generic mixed states. In Sec. [3] we 
consider the high temperature expansion of the partition function and compute the first 
moments of the purity, by making use of Zuber's solution of some basic integrals over 
the unitary group [TTJ. In Sec. [4] we establish the connection between our problem and 
the properties of the twirling maps. We conclude in Sec. [5] by summarizing our findings 
and discussing them in terms of future perspectives. 

2. The statistical approach: partition function 

In this section we review the statistical approach introduced in Ref. [9] for studying the 
local mixedness properties of pure states and discuss its generalization to the case of 
mixed states. 

Consider a bipartite system X = AB made up of two subsystems A and B, 
described by the Hilbert space Tix = Ha <8> W-b, with dim 7^ = Na, dim Kb = Nb 
and dimHx = N = N^Ng. Without loss of generality we will assume that Na < iVg. 
The states of X are represented by the set (5("Hx) of nonnegative unit-trace operators 
(density matrices) on the Hilbert space %x- The purity of such states, defined as 

7r AB (p):=Trp 2 G[l/iV,l], (1) 

for each p G &(1-Lx), provides a characterization of the global mixing of the system 
and induces a partition of G{%x) into a collection of distinct subsets & x (Hx) '■= {p £ 
&(T-Lx) '■ Trp 2 = x}. The minimum value of x — 1/N is attained when X is in the 
completely mixed state I/N, whereas the maximum x = 1 is attained over the set 
6i(Kx) consisting of all pure states \ip)x- For each p e <S("Hx) we can also define its 
A-local and B-local purity functions as 

7r A (p) := T^ , 7i B (p) := Trp| , (2) 

with pa = Tr B p and p# = Tr^p being the reduced density matrices of the subsystem 
A and B, respectively. 
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2.1. Total system in a pure state 

On the special set &i(T-Lx) of pure states p = l^x^l °f X, the A and the 5-local 
purities coincide 

k{1>):='ka(\1>)x(1>\)='Kb(\'iJ>)x(''!>\), (3) 

and provide a measure of the bipartite entanglement between A and B: the smaller 
7r(^), the larger the entanglement contained in \if))x(?P\- The statistical distribution of 
on &i(Hx) has been studied in Refs. [HI ED]. This was done by introducing the 
partition function 

Z(fi)= [ d^We-^W, (4) 



where the local purity tt(i[}) of \ip)x plays the role of an effective energy of the system, f} 
is a Lagrange multiplier that fixes the value of the purity/energy and selects an isopurity 
manifold [T2], and dpijp) is a (normalized) measure on the space of pure states &i(T-Lx)- 
The natural choice for the latter is induced by the Haar (probability) measure dpn(U) 
on the unitary group Ui^ix) — U{N), through the mapping 

\^)x:=U x \^o)x , (5) 

with \ipo)x an arbitrary reference unit vector of Hx- Thus the partition function 
becomes 

Z{p) = f dp H (U x )exp (-/3Tr(Tr B (f/ x ^ )x^ol4) 2 )) • (6) 

Then, by noticing that for every Ux £ U{N) the reduced density matrix p A = 
Ttb(\iP)x (01) can be written as pa = U A A A U A , with Ua £ U{1S 'a) and = 
diag{AA,i, ^a,2 ? • • • ^a,n a }, the expression ^ becomes 

Z{fi) = J dp H (U A ) J da(A A ) e^ TrA - = J da(A A ) e~^ A \ (7) 

where we exploited the fact that the measure induced over the density matrices pa 
by dpuiUx) factorizes into the product of a measure over the unitary group dpH(U A ) 
(related to the eigenvectors of p A ) and a measure dcr(A^) over the (N A — l)-dimensional 
simplex of its eigenvalues X A j [13J, |TJ] . In particular, it can be shown that [EJ [T71 [18] 

da(A A ) = C Na , Nb 8 ( 1 - A ^J II d ^)^AT NA 

\ l<i<N A J l<j<N A 

X J] (^A,l-XA,mfd NA \A, (8) 
l<l<m<N A 

with 

Na ' Nb rWx-i nN B - j)r(N A -j + iy w 
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and 8(x) and T being the unit step and the Euler gamma function, respectively. The 
derivatives of 2((3), evaluated for (3 = 0, yield the moments of ir(ip) with respect to the 
measure dfi(ip), i.e. 



((T^De^) = / <W) irW" = (-1)" ^ 



■ (10) 

p=o 

These functions fully determine the statistical distribution of 7r(i/}) on &i(Hx) and, in 
the high temperature regime, provide an expansion of Z((3). More generally, in analogy 
to what is commonly done in statistical physics, one can also define the moments of 
7r(V>) for (3 7^ as 

MM := ((Trpl) n )^ i(nx) = J d^) = ^ , (11) 

with d/j,/3(ip) being the canonical measure 

e -/Mv) 

dfipty) := d/i(^) -^r- . (12) 

The latter is a deformation of the Haar measure dfi(ip) obtained by including a 
nonuniform weight which explicitly depends upon the local purity, through j3. In 
particular, as (3 increases d/ig (■?/>) enhances the role of the states with lower values of 
Tr(ip) (i.e. larger values of bipartite entanglement), to the extent that for — > +oo only 



the maximally entangled elements of &\{Hx) contribute to the values (11). Since 7r(ip) 
is bounded above, one can also consider negative temperatures, where the role of less 
entangled states is enhanced, and when (3 — > — oo only the separable elements of &i(%x) 



contributes to the values (11). Consequently, across different ranges of temperatures, 
the moments in Eq. (11) characterize the statistical distribution of the local purity 
of X computed with respect to a canonical ensemble whose constituents are selected 
according to an effective thermal distribution characterized by the parameter (3. 

In the limit of large N, the /3-dependence of the statistics of the local purity n(ip) 
(and hence of the bipartite entanglement of the system) was characterized in Ref. [U [10] 



by identifying the class of states which maximize the distributions (12), i.e. typical 



states with respect to the canonical measure (12). In this context it was shown that the 
system undergoes two main phase transitions, related to different distributions of the 
eigenvalues of the typical states: a second-order phase transition, mentioned above, 
associated to a Z 2 symmetry breaking, and related to the vanishing of some eigenvalues 
of pa, followed by a first-order phase transition, associated to the evaporation of the 
largest eigenvalue from the sea of the others. 

2.2. Total system in a mixed state 

A natural question is what happens when the global system X is in a mixed state p of 
purity x < 1, rather than in a pure state. A generalization of Eq. Q is obtained by 
replacing z(ip) with (say) the A-local purity vta(p) of Eq. ^ and the measure d/j,(ip) 
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with a proper measure dp x (p) on the set & x (Jix)- This yields the following definition 
of the partition function of the A-local purity 

Z A (x,P) = [ dp x (p) e -P^[p) =Cx f dyU ( p) 5(Trp 2 _ x) e -pn A ( P ) ) (13) 



where dp(p) is a probability measure on the set of mixed states (see below), and 
C x '■= [J d/i(p) 5(Tr p 2 — x)] 1 a normalization factor. An analogous expression for the 
5-local purity partition function Z B (x, (3) is obtained by replacing 71a(p) with ttb(p) in 



Eq. (13): notice, however, that, at variance with the case analyzed in Ref. (9j HO], for 
x < 1 the partition function Zb{x, 0) will in general differ form Za{x, 0). 

It is worth stressing that the function Za{x, (3) provides only statistical information 
on the local mixedness of X, but not directly on its bipartite entanglement properties: 
this is due to the fact that for generic mixed states p of X the local purities ha{p) 
and txb{p) are not entanglement measures. A generalization of Eq. Q that retains 
the ability of characterizing the statistical properties of the bipartite entanglement of 
X for x < 1 could in principle be constructed by replacing ir(ip) with the convex- roof 
counterpart of tt a (p), namely 7Ta(p) = max^ J2jPj ^A^j), where the maximum is taken 
over all ensembles £ := {pj, \ipj)x}j which yield a convex decomposition of the mixed 
state p (see e.g. Refs. [201 E2])- The quantity 7Ta(p) is a proper measure of the bipartite 
entanglement, but the resulting partition function does not allow for a simple analytic 
treatment and will not be discussed in the present paper. 

Finally, since there is no unique measure on mixed states [15] ; we need to properly 



specify the choice of dp(p) which enters in Eq. (13) - the consistency requirement that 



for pure states Eq. (13) should reduce back to Eq. Q (i.e. Z A (1, (3) = Z((3)) does not 
eliminate such ambiguity. Indeed, as previously discussed for pA, the Hermitian matrix 
p can always be diagonalized by a unitary operation, and as consequence we can write 
the measure dp(p) as the product of a measure on the (N — l)-dimensional symplex of 
the eigenvalues and a measure on the unitary group U(N) related to the eigenvectors. 
However, if on the one hand it would be natural to take the Haar measure on U(N) 
so that d/i(p) = dpjQJxpUx), on the other hand the measure on the eigenvalues can be 
chosen in different possible ways [23] . In order to overcome this ambiguity we will use 
a balanced purification strategy. Let us introduce hence the composite Hilbert space 
Hix' = Hx ® Hx', where %x — Hx 1 are isomorphic. In this iV 2 -dimensional Hilbert 
space, each p of X can be represented by those pure states \^)xx' which provide a 
purification for such density matrix, i.e. which satisfy the identity 

p = Tr x ,(\y} xx ,(V\) . (14) 

Thanks to this identification we can now induce a measure on (3("Hx) by sampling the 
pure states on Hxx' according to the unique, unitarily invariant Haar measure which, 
as usual, is induced by the Haar measure on the unitary group U(N 2 ) through the 
mapping \^)xx' '■= Uxx>\^o)xx', where \^o)xx' is an arbitrary reference vector and 
Uxx> & U(N 2 ). With this choice the partition function becomes 

Z A (x,/3) = C X J dp H (U xxl ) 5 (x - Tr(Tr x mxx'(n 2 ) e~ fi H^*'^™ 2 ), 



Statistical distribution of the local purity in a large quantum system 



7 



(15) 

where we used the fact that p A = Tr B p = Tr B (Tr X i\^)xx'(^\)- Analogously to 
what we have seen for the pure case, x = 1, by writing p = UxAxU x w hh Ax = 
diag(Ai, A 2 , X N ), we get 

Z A (x,P)=C x j dp H (U x ) J da(Ax) s(x-TrA x ) e -P-*{V*<P*W*W) , (16) 

where dpn(Ux) is the Haar measure on 14 (N) and 

da(Ax) = c N 5 (i - a n n ( a 

V l<i<N J l<i<N l<i<j<N 



\ 3 yd N \, (17) 

<3t 



with 



r r(iv 



r(iv + i)ni< fc <ivW 



Therefore, we have identified the measure dp x (p) of Eq. (13) with 



dp x (p) = C x dp H (U x )da(A x )6(x - Tr A^) . (19) 

Notice that in the case of pure states, i.e. x — 1, the density operator of the system 
reduces to p = \ijj)x($\, where \ip)x = Ux\?Po)x, \ipo)x being an arbitrary reference state 
(see Eq. (|5p), and the matrix Ax becomes a rank one projection. Thus, the expression 
(16) reduces to ([6]), namely 

Z A (\,p) = Z{p). (20) 
2.3. Asymptotic behavior and analysis of moments 



For x < 1, the integration over the unitary group 14 (N) in Eq. (16) does not factorize, 
making the computation of the partition function far more complicated than for the 
case of a pure state ([7]). The only notable exception is the case of maximally mixed 



states (i.e. x = 1/N), when the Dirac delta in Eq. (16) selects a unique diagonal matrix 
Ax (the totally mixed state of X). This makes the exponent equal to e _/3//7VA for all Ux 
and yields the following exact expression 

Z A {l/N,0)=e-f" N *. (21) 

Otherwise, for intermediate values of the purity, 1/N < x < 1, the situation is much 
more complicated. Still, as we will show in the following, at small the evaluation of 
the moments A4^(x, (3) associated with Z A {x, (3) admits an exact analytical treatment. 
The latter are formally defined as 

Mi( X , 0) := / W ) , 5 W = jglj ™0® , (22) 

and represent the average value of 7r^(p) with the canonical measure 

g-/3 vr A (p) 

dp x Ap) ■= &Vx(p) ^7 — wr , (23) 
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with dp x (p) given by Eq. (19). For pure states (x = 1) the Ai^(x,/3) coincide with 



the moments Ai n (/3) defined in Eq. (11): at f3 = 0, in the large N limit the expression 
for such quantities has been computed in Ref. [9], while the exact expressions for first 
five of them can be found in Ref. |24j. In the case of a totally mixed state (x = 1/N) 



Eq. (21) yields instead values which are independent of the temperature (3, namely 

Mi{l/N, /?) = Mi{l/N, 0) = N A ~ n . (24) 



For intermediate values of x, by expanding Eq. (22) up to the first order in (3, we get 



M£(x, 0) ~ M*(x, 0) - p [M£ +1 (x, 0) - Mf(x, 0) M*(x, 0)] , -> 0. (25) 



Incidentally, notice that in agreement with Eq. (|24j), the /3-corrections of Eq. (25) vanish 
when x = 1/N. The above expression shows that, at least in the high temperature 
regime, we can focus on the unbiased moments Ai^(x, 0). 

3. Moments of the purity at f3 = 

3.1. First moment 

In this section we compute the exact first moment of the purity Aif(x, 0), by making use 
of Zuber's solution of some basic integrals over the unitary group [11]. In particular, we 
will show that the only dependence on the spectrum of the density matrix of the global 



system is in terms of its purity x, whose value is fixed in the partition function (13). 
Let us fix the spectrum of the global density matrix p of the system: 

A x = diag(A Ars(Q _i )+/ g), (26) 

with a = 1, . . . , Na and (3 = 1,..., Nb- A purification of Ax in the space Tixx 1 , with 
X = AB and X' = A'B' is 

\*)xx' = J2 J2 WP)ab ® \aP) A > B >, (27) 

l<a<N A 1<P<N B 

where we have set, for simplicity, X al3 = X Nb ^i^ +/3 and \ci/3)ab — \°)a ® \0)b, {1°)^} 
and {\(3)b} ({\ a )A'} and {\(3)b'}) being the reference basis in Ha and Hb (Ha 1 and 
Hb 1 ), respectively. The set of vectors in Hxx' with the same Schmidt coefficients is 
given by Uxx'\^)xx', where Uxx 1 = Ux <8> Ux>, with Ux,Ux> E U(N), and yields 
the set of density matrices with the same spectrum Ax, namely p = UxAxU x . By 
partial tracing over subsystem B one obtains the set of reduced density matrices 
Pa{U) = Ttb(^x'(Uxx'\^)xx'{^\U xx ,)). Notice that this expression does not depend 
on U x > eU(N): 

Pa= ^i3^B(U A B\a(3) A B(a(3\U AB ) 

l<a<N A 1</3<N B 

= J2 X ^B(j\U A B\af3)AB(af3\U A B\j)B , (28) 

1<q<JV a l<P,j<N B 
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k A {U x A x U x ) = Trp 2 A = Yl Yl A «i/3i^ 2 /3 2 

l<ai,a 2 <N A 1</3i,/3 2 <Nb l<h,h<N B 
X AB («2/32 1 U\ B \j 2 ) B ( jl I U A B I OtiPi ) AB 

X AB{0ilPl\U\ B \ji) B {j 2 \U AB \0i2P2)AB , (29) 

which, by the completeness relation for subsystem A, becomes 

ir A (U x AxU x ) = Y Y Y Y •VA'Wa 

l< ai ,a 2 <N A l</3i,0 2 <N B l<h,i 2 <N A l<ji,j 2 <N B 

X AB^^lUAB^lh) AB(iljl\U A B\aiPl) AB 

X A B{0Hpl\U AB \i 2 ji) AB {i2j2\U A B\0l2p2)AB- (30) 



Let us now compute the first moment of the purity (22) at f3 = 0. By recalling 
that Z A (x,0) = 1, dp, xfi (p) = dp x (p), and dp x (p) = C x dp H (U x )da(A x )d{x - Tr A^J, 
we get 



where 



M£{X,0)= j dp x ( P ) 7T» 

= C X j da(A x )S(x - Tr A 2 X ) M^(Ax) =: (M^(Ax)) x , (31) 
Mt{A x ) := j dfi H (U x ) ix n A {U x A x U x ). (32) 



From (30) the average over the unitary group of the first moment particularizes to 
Mi(A x ) — ) j A ai(gi A Q2/ g 2 

l<ai,a 2 <N A l<Pi,P2<N B 

x Y Y / d ^(^)^iii.«i/ 3 i^2i2,«2/3 2 ^i 2 /3 2 ,uj 2 ^ 1 /3 li i 2 j 1 ) (33) 

l<h,i2<N A l<ji,j 2 <N B J 

where Uij >a p = ab^jWabW^ab- This integral can be explicitly computed by using 
Zuber's solution [TTj : 

I dp H {U)U lin ...U ln]n {U)Ut lh ...Ul ln = CM II ^ a ,£ T{a) )5(j a ,k Ma) ),(3A) 



T,a£S n l<a<n 



with 



n , (x (fc) (i))V fc) (M) ^ 

|y|=n ftV y 

where C[a] is the sum over the Young diagrams Y of the character x (M) °f the 
symmetric group associated to Y, depending on the conjugacy class [a] of the 
permutation a, Sfc(J) is the dimension of the representation, and 5 (a, b) is the Kronecker 
delta. Applying this solution to ( J33| ) we get 

My(A x ) — ) j A ai/3l A Q2ft , 

l<ai,a 2 <A r A l<Pi,P2<N B 
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X Y C '[ "]/l( r )^( a l^l' a ^(2)/ 3 T CT (2))5(a2/32,ar CT (l)/3T CT (l)) 

Y Y Y Y c, M/i( r ) 5 ([ T(7s ]' c ) A «c(i)^(i) A « e(2) ^ (2) , 

l<ai,a 2 <N A l<Pi,fa<N B T,a&S 2 cGC(5 2 ) 

(36) 

where /i(vr) depends on the permutation n G S 2 

M 71 ") : = Y <Kn>Ml))<K«2,M2)) Y S (h,j7v(2))S(j 2 , (37) 
l<ii,i 2 <^ l<ji,j2<N B 

s G S*2 is the transposition (swapping) of pairs of nearby indices ([s] = [2]) 

= «2 and i s{2) = h , (38) 
and C(S 2 ) = {[l 2 ], [2]} is the set of the conjugacy classes of the symmetric group S 2 . 



From (36) it can be easily inferred that the only possible contributions of the spectrum 



are related to the conjugacy classes of the symmetric group S 2 : 



[ras] = [I 2 ] 



rers 



[2] 



Y Y Xa ? 

yl<a<N A 1</3<N B 

E E ^ 

^l<a<N A 1</3<N B 



TrA^. 



(39) 



By summing and by using the explicit expressions of the coefficients ( [35] ) [TT] 

C[! 2 ] = ,» r ,Lr ■ TT ' C[2] = - — ■ (40) 



(JV-l)(iV+l)' 



(AT — 1)N(N + 1) ' 



we get 



M A fA AT A (iV 2 -l) iV^iV 2 -!) 



(41) 



The first moment of the purity of subsystem A, is the average of (41) over the spectrum 



of the system. By plugging fl41ft into (31), we finally get 



Mf(x,0) 



N A (N 2 B -1) t N B (N 2 A -1 



+ 



N\Nl - 1 N\Nl - 1 



-x. 



(42) 



Notice that for x — 1 this expression reduces to the one given in [15], [16] (TT] [191 EH EH] . 
Some special cases:- It is worth noticing that for a balanced bipartition = Nb = 



N > 1 Eq. (42) yields 



Mf(x,0) 



N(l + x) l + x 



N+l 



N 



(43) 



At x = 1 (i.e. pure global states), Eq. (43) coincides with that obtained in 



Refs. [HI [E3 El]- Finally, consider the case in which p is maximally mixed, i.e. is 

(44) 



the density matrix I/N. In this 1/N and Eq. (|42[) gives 



M?(1/N,0) 



N A 



in agreement with the general result (24). 
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3.2. k-th moment 

The technique shown in the previous section can be easily generalized in order to 



compute from (|32|) higher moments at (3 = 0. We get 
dfi H {U x ) n k A (U x A x U x ) 



Mf(A x 



l< ai ,...,a 2k <N A l<Px,...,fc k <N B l<i<2k 



l<h,-,i2h<N A l<ji,-,hk<N B 
t rrt 



e / d ^(^) n ^ 



!<£<2k 



X 



nTTt r/l 
m — 1 )^2mj2m- 1 



(45) 



Km<k 



Equation (34) for n = 2k gives 

-m^(Ax)= £ E E E w*o-w[tH, c ) n A 



a c (i)Pc(j) ' 



l<an,...,a 2fc <JVA l</3i,...,/3 2 fe<A r s T,<re5 2 fe ceC(S 2fe ) 



Ki<2fc 



(46) 



where /fc(7r) depends on the permutation ir e 

fk(n):= ^ Ell <K^'M<o) JI 5 0'2m-l,J 7 r(2m))^(j2m, j7r(2m-l)J 



l<il,...,i 2 *<JVA l<jl,-,hk<N B l<£<2k 



Km<k 



and, analogously to Eq. (38), s is the swapping of pairs of nearby indices 
i a (2£-i) = ^ and 



(47) 



(4£ 



« a (20 = ke-i V£ = l,...k. 
Observe that when k = 1 we retrieve M.f(A x ) (see Eq. (36)). The different 



contributions of the spectrum can be classified in terms of the the conjugacy classes 



of the symmetric group, as shown in Eq. (39). However, for k > 1, they do not 
depend only upon the purity x = TrA 2 ^, but exhibit a more complex dependence on 
the spectrum, through its higher order invariants Tr A^, with k > 2. Thus the integral 



on the spectrum (31) is in general non trivial. 



3.3. Second moment 

Now let us fully compute the second moment of the purity for arbitrary bipartite 
states, with purity x G [1/N, 1], generalizing some results found for pure states, x = 1, 



[HI [151 El]- The second moment can be directly computed by setting k = 2 in Eq. (45) 



The expression for the coefficients C[ir] in (34), when tt 6 C(Si) is [TT] 
c [ l4 ] - 8iV 2 + 6 

C [2, l 2 ] 
C[2 2 ] 



(N - 3)(JV - 2){N - l)N 2 (N + 1)(JV + 2)(iV + 3) 
1 



(JV - 3)(iV - l)JV(iV + 1)(JV + 3) 
N 2 + 6 

(N -3)(N- 2)(N - 1)N 2 (N + 1)(N + 2)(N + 3) 
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(N - 3)(JV - 2)(JV - 1)A^ 2 (A^ + 1)(N + 2)(N + 3) 

5 



(N-3)(N-2)(N-l)N(N+l)(N + 2)(N + 3)' v ; 
The symmetric group S 4 , consists of five conjugacy classes, that yield the following 



contributions to the integral (46) in terms of the spectrum of p: 



\Tcrs\ 



[1 4 ] 



E E Xc *p 

a<a<N A 1</3<N B 



[t<ts] = [2, l 2 ] 



H = [2 2 ] 



[rcrs] 
[ras] 



[3,1] 
[4] 



E E 

vl<oi<AfA l<Pi<N B 

E E 

v l<o<AT A 1</3<7V S 

E E 

v l<oi<Af A l<Pi<N B 

E E 

l<a<Af A 1<0<N B 



X 2 




X 3 



( E E 

\l<a 2 <N A 1</3 2 <N B 

(TrA?.-) 2 

( E E 

\l<a 2 <N A 1</3 2 <N B 



A 



"2/92 



A 



<*2/32 



A a/3 



TrA 



with r,(T 6 54 and s £ S 2 defined in (|48|) being k = 2. By gathering all we get 



C[1 4 ]A^V 4 + C[2, 1 2 }2N A (N 2 A + 2) N% 

+ C[2 2 }N 2 A (N 2 A + 2)N 2 B + C[3, l]8iV 2 iV| 
+ C[4]2iV A (2iV 2 + 1) jV fl 

+ Tr A 2 X \C[1 4 ]2N A (N 2 A + 2) N 3 B 

+ C[2, 1 2 ]2N 2 A N 2 (N 2 A + 3N 2 B + 14) 
+ C[2 2 ]2N A N B (N 2 (N 2 + 4) + 2 (N 2 + l)) 
+ C[3, l]8jV A jV B (iV 2 (iVl + 2) + 2iVl + 1) 
+ C[4]4A^ 2 (N 2 A + 8) iV| 

C<[l 4 ]iV 2 (N 2 A + 2) iV 2 

+ C[2, 1 2 }2N A N B (N 2 A (N 2 B + 4) + 2 (iVl + 1)) 
+ C[2 2 ]N 2 A N 2 (2Nl + 3N 2 + 4) 
+ C[3, l]24A^i 

+ C[4]2iV A iVB (2JVj[ (iV| + 1) + 4A^| + l) 



(TrA|) 2 



+ TrA 



x 



C[1 4 ]8AT 2 ^ 



B 



+ C[2, 1 2 ]8N A N B (N 2 A (N 2 + 2) + 2N 2 B + l) 
+ C[2 2 ]24A^Vf 



TrA 2 , 



TrA^ 



(50) 
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+ C[3, l]SN 2 A N 2 B {N\ + N% + 6) 

+ C[A]8N A N B {N\ (N 2 B + 2) + 2N 2 B + l) 



13 



+ Tr A 



x 



C[l 4 ]2N A (2N 2 A + 1) N B + C[2, l 2 ]4iV 2 (]V 2 + 8) N 

+ C[2 2 ]2N A N B (2N 2 A (N 2 B + l) + AN 2 B + l) 
+ C[3, l]8iV A iV B (iV 2 (JV| + 2) + 2iVl + l) 

+ C[A]2N 2 A N 2 B (N 2 A + 3N 2 B + 14) 



(51) 



from which it follows 

Mi(k x ) = c NA)N 



(N 2 B - 1)(N A N 2 B (N 2 B - 1) - 2N 2 A (6N 2 B - 7) + 22) 

+ Tr (2N A N B (N 2 A - 1)(N 2 B - l)(N 2 A N 2 B - 14)) 

+ (Tr A^) 2 {N 2 A - 1){N B N A + N B N A - UN 2 A N 2 B + 6N 2 B + 30) 

+ Tr A 3 X A0(Nl - 1)(JV| - 1) 

+ TrA^ (-10N A N B )(N 2 A -1)(N 2 - 1) 



where 



1 



CN A ,N E 



(52) 



(53) 



NlNl{N\Nl - 7) 2 - 36" 

This expression generalizes the already known result for the pure case, when Tr A^ = 1 
for all k [IU US]. In particular if N B = N A = y/N we get 

M^(Ax) = | [C[1 4 ]AT 3 + C[2, 1 2 ]2N 2 (N + 2) 

+ C[2 2 ]iV 2 (A^ + 2) + C[3, 1]8N 2 + C[4]2N(2N + 1) 



+ Tr A 



A' 



C[1 4 ]2AT 2 (N + 2) + C[2, 1 2 ]AN 2 (2N + 7) 
+ C[2 2 ]2N(N(N + 6) + 2) 
+ C[3, 1]8N(N(N + 4) + 1) + C[A]AN 2 (N + 8) 

+ (Tr A^) 2 [c[l 4 ]AT 2 (Af + 2) + C[2, 1 2 ]2N(N(N + 6) + 2) 

+ C[2 2 ]N 2 (5N + 4) + C[3, l]24iV 2 + C[A]2N(2N(N + 3) + 1) 

+ Tr A^ [c[l 4 ]8N 2 + C[2, 1 2 }8N(N(N + 4) + 1) 

+ C[2 2 ]2AN 2 + C[3, 1}1QN 2 (N + 3) + C[A]8N(N(N + 4) + 1) 

+ Tr A^ \c[l 4 ]2N{2N + 1) + C[2, 1 2 ]AN 2 {N + 8) 

+ C[2 2 ]2N(2N(N + 3) + 1) + C[3, 1}8N(N(N + 4) + 1) 



C[4]4iY 2 (2N + 7) 



(54) 



with 



M*(A X ) = c N 
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that is 

(N 5 - 2N 4 - 11N 3 + 2QN 2 + 8N - 22) 

+ Tr A 2 X (2N 5 - AN 4 - 26N 3 + 5QN 2 - 28N) 
+ (Tr A^) 2 (N 5 - 15N 3 + 20N 2 + 2AN - 30) 
+ TrA^ 40(iV- l) 2 

+ TrA^ (-IW)(N - l) 2 
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cn 



N 2 (N 2 -If -36 



(55) 
(56) 



In the thermodynamical limit, N » 1, we find 



Mi 



N 



{l + x) 2 + 



1 

AT2 



(57) 



From Eqs. (42), (52) and (53) we can now compute the exact expression for the 
second cumulant of the purity at (3 = 0: 



^(x,0) = M£(x,0) - (M?(x,0)Y 



+ 



2(N 2 A -1)(N 2 



(N 2 A N 2 - 1) 2 {N\N% - 13N 2 A N 2 + 36) 



+ x 



2(N 2 A - 1)(N 2 - 1)(-2N A N B )(N 2 A N 2 + 11) 



+ x 



(N 2 A N 2 - 1) 2 {N\N% - \3N\N 2 B + 36) 
2 2(7V 2 - l)(iV 2 - 1)(N A N B - 4N 2 A N 2 + 15) 



(N 2 a N 2 b 



l) 2 (N\N B 



13N 2 A N 2 B + 36) 



+ (TrA^ 

+ <TrA^ 
See Ref. [151 121] for the case x = 



A0(N 2 -1)(N 2 



1) 



N 2 A N 2 (N 2 A N 2 - 7f - 36 
(-10N A N B )(N 2 A -1)(N 2 -1) 



N 2 N 2 {N 2 A N 2 



A ly BK ly A ly B T) 2 

when all the traces are 1. 



36 



(5* 



3.4- High temperature expansion of the first moment of the purity 

We can now compute the approximate expression for the first moment of the purity for 
small f3, by plugging Eqs. (Elk and (58) in Eq. (25), with n = 2: 



M?(x, (3) ~ M?(x, 0) + [(M?(x, 0)) 2 - M$(x, 0)] 



N B (N 2 A -1) N A (N 2 -1) 



N\N 2 B - 1 



N\N 2 B 



+ 



2(N 2 A -1)(N 2 



1)(N 2 A N 2 + 11) 



(N 2 A N 2 - 1) 2 {N\N B - \3N\N 2 B + 36) 



x 



2{N\ 



1)(7V 2 - l)(-2N A N B )(N 2 N 2 + 11) 



(N 2 N 



B 



lf{N\N% 



13N\N 2 + 36) 
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Figure 1. (Color online) First moment of the purity as a function of 0. The horizontal 
(red) line refers to the set of totally mixed states, with Mi(l/N,f3) = 1/Na; the 
(black) curve refers to pure states for (3 > 0; the vertical (blue) line corresponds to 
the first moment of the local purity at /3 = 0, for arbitrary mixed states; finally, the 
shaded (light-blue) region refers to high temperatures. 



X 



2{N\ - lXjYj - 1){N\N% - AN 2 N 2 + 15) 
(N 2 A N 2 - l) 2 (iV 4 iV 4 - UN\Nl + 36) 



m{N\-i){Nl-\) 

N\Nl{N\Nl - 7) 2 - 36 
{-IW A N B )(N 2 A -1){N 2 



N\Nl{N\Nl - I) 2 - 36 
For a balanced bipartition, N B = N A = yN, we get 

M?(x, 0) ~ A<tf (a, O) + [{M?{x, 0)) 2 - jVtf (a;, 0)] 



x 



1 - 
+ P 



N 



2(N 2 + 11) 



+ x 



(N + l) 2 (iV 4 -13Ar 2 + 36) 
AN(N 2 + 11) 



— x 



(N + l) 2 (N 4 -13N 2 + 36) 
2 2(iV 4 - AN 2 + 15) 
(iV + 1) 2 (N 4 - 13N 2 + 36) 
40(N- l) 2 



-(TrA^) x 
- (TrAi), 



N 2 (N 2 - 7) 2 
(-10JV)(JV- 
N 2 (N 2 - 7) 2 - 36 



-36 
D 2 



and in the thermodynamical limit 

1 + x 



M?(x,f3) 



2/3 



N N 2 



x 



o 



AT3/2 



(59) 



(60) 



(61) 
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One might wonder whether higher order cumulants follow a pattern similar to (61) 



Notice that (61) suggests a convergence radius for the high temperature expansion 
fi c ~ iV 3 / 2 (l+x)/2x 2 , which grows indefinitely when x — > 0, in accordance with Eq. ( |2l[ ). 
See also Fig. [T] 



Equation (61) at x = 1 can be compared with the results of Ref. [D] where f3 was 



replaced by the scaled quantity (3 = (3'N 3 ^ 2 . With this choice our expression yields 



M?(l, /3'iV 3 / 2 ) ~ Mf(l, 0) + P'N 3 / 2 [(M?(l, 0)) 2 



N 



(62) 



in perfect agreement with the behavior reported in Fig. 2 of Ref. [9] . Figure [T] yields an 
overview and summarizes our results. 



4. An alternative approach based on Twirling 

In this section we will establish an interesting connection between our problem and 
the theory of quantum channels. In particular, we present an alternative approach for 



computing the moments A4^(x,0) of Eq. (22) which exploits the properties of twirling 
transformations |26l 1271 1281 [201 [21] . 



For explanatory purposes we start in Sec. |4.1| by deriving a general expressions 
for the moments M^(x = 1, 0) associated with the case in which p is pure (i.e. the 



quantities M. n of Eq. (10)) and verify that it yields the exact value given by Lubkin [15 
for n — 1. The case of mixed states is then addressed in Sec. 4.2 showing that for n = ] 



it reproduces the results of Sec. 3.1 



4-1. Pure initial states 

Let us consider a fixed (normalized) pure state of the global system X = AB, 
\i^o)x, and parametrize the pure states of Hx as in (|5f, i.e. \ip)x '■= Ux\ipo)x, with 
Ux € U(Hx) — U{N) distributed according to the Haar measure dpuiU). Its local 
purity ([3]) can be expressed as 

7T(V) = Tr (B(AU A B\A)AB(A\Ui B \t)B B^\UAB\^)AB{MU\B^)B y 

1<£,£'<N B 

= AB ^\ ( U ab\Z)b(?\U ab ) \^) ab ab (M {y{ B \i') B {l\U AB ) |^ )ab , (63) 

i<e,i'<N B 

where is an orthonormal basis ofH B , and the cyclicity of the trace was used. We 

can recast this expression into a more compact form by doubling the Hilbert space, i.e. 
adding two auxiliary copies A' and B' of A and B, respectively. We get 



7r (■?/>) = Tr 



(Uab ® Ua>b>) ( \^o)ab(^o\ <8> |V'o)a'B'(V , o| 
x (S B \b> ® Iaa') 



U 



t 

AB 



A'B' 



(64) 
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where the trace is over all degree of freedom (i.e. AA'BB'), Iaa 1 is the identity operator 
on A A' and 



s B \b>-= \t)B(e'\®\e')B>(e\, (65) 

i<e,e><N B 



is the SWAP operator on BB' - this is the unitary, self-adjont transformation which, 
for all operators 0# and T B >, gives 

S b \b>{®b®Tb')Sb\b> = T b ®Q b ,. (66) 

Remembering that first moment Mf(x = 1,0) = Mi of Eq. (10) is obtained by 
averaging over all possible U x , we can then write 

Mi = Tr (r (2) (|^o)x(^o| ® \i>o)x'(H) {Sb\ b> ® Iaa>) , (67) 

where X = AB, X' = A'B', and where is the Completely Positive, Trace Preserving 
(CPTP) twirling channel [261 EH EH E2] which transforms the operators ®xx* of XX' 
into 

T {2 \Qxx>) = I dfi H (U) (U x ® U x .) G X x> (U x ® (68) 



This map plays an important role in quantum information theory where it was first 
introduced as a tool for characterizing the distillability of bipartite entanglement [2TT 120] . 
It has several properties which allows us to simplify the calculation. For instance 
it is known that maps all the states of the system into (generalized) Werner 
states (2ni EZ]- Furthermore it is self-adjont - i.e. its description in Heisenberg picture 
coincides with T*- 2 **. In particular this last property can be used to rewrite (67) as 

Mi = (x{M ® x>(H) T {2) (S B \b> ® Iaa>) (l^o)x ® \ih>)x>) ■ (69) 

Explicit expressions for the action of can be obtained by exploiting the symmetry 
of dfj,H(U). In particular it is possible to show that T^ 2 '(Q X X') can be decomposed as 
a linear combination of the projections on the symmetric and anti-symmetric subspaces 
of XX' = ABA'B' (with respect to the bipartition AB\A'B'). Introducing then the 
SWAP operator which exchanges X with X' this can then be written as 

t-(2)/ \ NI X x' - S x \x' „ „ NS x \x> - Ixx> „ x , . 

T {QxX,) = iV(iV2 _ i) Tr Qxx ' + iv(iV2 _ i) ^( S x\x>®xxi) (70) 

NTrQ XX i - Tr(S X \x'9xx') f NTr(S X \x'9xx') ~ Tr9 X x' „ 

N(N 2 - 1) XX ' N(N 2 - 1) XlX ' ' 

(71) 

where = NaNb is the dimension of "Has- (Here Ixx' is the identity operator while 
Sx\x> = •S'abia'B' = Sb\b' ® iS^ia' is the SWAP which exchanges AB with A'B'). 

Thanks to this expression we can now easily compute the value of Mi, either 
using Eq. (67) or Eq. (69). Consider for instance the first approach. We have first 
to compute the quantities Tr Qxx' and Tr(Sx\x'®xx') with Oxx> being the operator 
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|V'o)x(V'o| ® |V'o)x'('0o|- This is 

Trd^x^ol ® iVoMVtol) = 1 , 

Tr ((S B \ B > <8> S'aiaOIV'oJab^oI ® l^oW(^ol) = 1 , (72) 
where in the second expression we used the fact that \ip) AB ® W)a'b' is invariant uder 
S x \x>, i.e. (Sb|b' ® SU[A')(I^)ab <8> I^A'S') = W)ab ® W)a'b>- Replacing all this in 



Eq. (70) we get 

T (2) (|^o>x^o|® l^o)x'(V'ol) 



Ixx 1 + S. 



X|X' 



N(N + 1) 



(73) 



and thus 



Mi = N ^ + Tr ((Iaba'bi + S B \ B > ® 5a]a0 (S'bib' ® 7u') 



N(N + 



(74) 



(here we exploited the fact that S B , B , = Ibb>)- Now we can use N = N A N B and the 
identities 



to get 



Tr (S B]B , <g> I AA ,) = N B N\ , Tr (l BB , <g> S A \ A >) = N A N B 
N A + N B 

Mi 



N A N B + 1 

which coincides with the correct value [T5J 121] ■ 

We mention that the same techniques can also be applied to higher moments M. 



(75) 
(76) 



The extension of Eq. (67) for n ^ 2 is obtained by introducing In copies of AB organized 



in the n pairs, i.e. AiB x A' x B' x , A 2 B2A' 2 B' 2 , ■ ■ •, A n B n A' n B' n . We then introduce the 
following generalized twirling transformation acting on XX' = A X B X A X B' X , A 2 B 2 A' 2 B' 2 , 
■ • ■, A n B n A' n B' n , i.e. 

T (2n) (©xxO = / dfi H (U) {U ® U <g> ■ • ■ ® U) ® xx , (W ® U ] <g> ■ ■ ■ <g> 17+), (77) 

J S v ' s v ' 

2n 2n 

with ©xx' being a generic operator on 7ixx> '■= 'H AB n - This channel is a proper 
generalization of the map whose properties can be established along the lines of 



r (2n)- 



Ref . [TTj . With this choice Eq. ( 10 ) then can be expressed as 

M*(x=l,0) = M n = Tr(V (2n) (|*® 2 ><^ 1 p.'s ■, , A 
where |\|/ (X,2 ) (g,n := ®™ =1 (\ip) AjBj ® \^)a'.b'^ Ia being the identity on the 2n copies of 



r 2 n }(S^ n ^ 



(7f 



3 j ' 

y(2re) 



A, i.e. .A = AjAi • • • A n A' n , and with Sg being the SWAP operator which exchanges 
B X B 2 ■■■B n with B[B' 2 ■ ■ ■ B' n pairwise, i.e. S§ n) = ®] =x S Bj]B ,. 

4-2. Mixed initial states 

Consider now the case with x < 1. Following the parameterization introduced in Sees. [2] 
and|3l we split the average over the set & X {T-L X ) of the density matrixes of global purity 
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x, as an average over the unitary rotations of acting on T-t x followed by an average 
over the space of the eigenvalues p - see Eq. (16). Specifically this is accomplished by 
writing p(U) = Uab Aab U ab , with Uab being a generic unitary transformation on 
Uab, while Aab represents a given arbitrary choice of the system spectrum, see Eq. 
(26). For convenience, let us rewrite the purification (27) of the density matrix p as 

\V)ABab= Yl Yl V 7 ^ |a>^® \P)b® \P) b , (79) 

l<a<N A 1</3<N B 

where now the ancillary systems, isomorphic to A and B, are labelled by a and b, 
respectively. The reduced density matrix pa(U) = Trg p(U) can thus be written as 

Pa{U) = Tr Bab{U AbY&) ABab 

J2 Bab{q\U A B\^)ABab{^\U ] AB \q) Bah} (80) 



B 



with {\q)Bab} being an orthonormal basis of Bab. Similarly the local A-purity of Pa{U) 
becomes 

7l A (p) = Trp 2 A (U) Tr (Bab(q\U AB \y)ABab(y\U AB \q)Bab 

l<q,q'<N A N% 

X Bab(q'\U A B\y)ABab(y\U AB \q)Bab 

= ABab(y\U AB \q)Bab(q'\U A B\y) ABab 

l<q,q'<N A N% 

X ABab{^\U ] AB \q) Bab{q\U A B\^) ABab (81) 

which, once more, can be casted as an expectation value on |\I/)® 2 by doubling the space 
(see Eq. (29)). Therefore by integrating over Uab we get, 

M?(A X ) = Tr (r (2) (|^)x,(^| ® |*>xv(*|) (S B ab\B'a>b> <S> 1 AA> ) ) 

= (xxi^l <8) X>xi(ytyT {2) (S B ab\B>aib> ® IaA'^) (\^)xx <S) \^}x'x'\ (82) 

where being the twirling transformation on XX' of Eq. (68) with X = AB, 
X' = A'B', x = ab, and x' = a'b' (here A',a',B',b' are the auxiliary copies of A,a,B, 
and b respectively). The above expression is the average purity of the subsystem A 
computed for states p having the same spectra A X - 

To compute the above quantity this time we use the last identity of Eq. (82). 
According to Eq. (71) we have to compute Tr Qxx' and Tr(Sx\x>®xx') with Qxx' 
being the operator S B ab\B'a'b' <8> Iaa* = S b \b' <8> S b \y <g> S a \ a ' <8> Iaa>- That is 

T^ABA'B'(SBab\B'a'b' ® I A A') = T^ABA'B' (S B \B' <S> S b \ b > <S> S a \ a > <g) I AA >) 

Tr(S B \B') Tt(Iaa') S b \ b > <8> S a \ a > = N 2 A N B S b \y ® S a \ a > , (83) 



TlABA'B' {SAB\A'B'{.S B ab\B'a'b' ® IaA')) ~ T^ABA'B' (IbB' ® Sa\A' ® S b \ b > Cg> S a \ a >) 

= Tt(Ibb') Tr(S AlA ') S bW ® S a{a , = N 2 B N A S b{b , ® S a{a , . (84) 
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Thus from Eq. ( 71 ) we get 

( S B ab\B'a'b' ® I AA' 



N B (N\-l) 

Ar2 Ar2 V~ BB' ® lAA' <£> db\V ® ^a\a' 

N\N B — 1 11 
Na{N 2 b - 1) 

N 2 N 2 — 1 6 I 6 ' a l a ' ' 



(85) 



where we used iV = N A N B . Replace now this into Eq. (82) and employ the identities 

Xx(V\ ®X'*(V\)(Ixx> ®Sx\x>)(\tyxx ® \y)x'x') = TrA x , 



The final result is thus 



Mf(A x ) 



N B (N 2 A -l) TrA ^ + N A (N B -l) 



N\N% - 1 



N\N% - 1 



57) 



which is Eq.(41) and depends upon the spectrum Ax only through its purity. By 



averaging upon on Ax while keeping fix x gives us the same result ( 42 ) . 



5. Conclusions 

From the results obtained in the previous sections one can infer that the same 
phenomenon of concentration of measure that occurs for the eigenvalues of the reduced 
density matrices of pure states [91 [10] occurs in the present case as well. Indeed, we 
observe that for large N the leading order of all fc-moments equals the fc-th power of 
the first moment (43), that, for a balanced bipartition Na = N B = v^V, reads 

:i+x) 



M?{x,0) 



N 



+ o 



AT3/2 



This observation spurs from the calculation of the second moment (55): 

{l+x) 2 



N 



+ o 



N 2 



Therefore, 



M${x,0) = Mf(x,0) 2 + O 



N 2 



(89) 



(90) 



By observing that the contributions of 0(1 /AT) to the second moment come from the 



coefficients proportional to C [l 4 ] in (46), i.e. the identity permutation a = 0, and 
conjecturing that the leading contribution for all the k-th moments comes only from 
C [l k ] ~ 1/N k one gets 

k 



Mi{x,0) 



1 



x 



N 







1 



AT(fc+2)/i 



(91) 



Another check of the validity of (91) derives by the interpolation between maximally 



mixed and pure global states. See Fig. [T} The scaling with A" is preserved, therefore 
allowing to interchange the x — > 1 and A" — > oo limits. 
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This uncovers the issue of computing the subdominant terms in the k-ih moments 
expansion that become instead leading order terms in the k-th cumulants. This could be 
possible if we had an appropriate asymptotic expansion of the combinatorics coefficients 
C, which we do not know. We leave this as a challenge for future work. 
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